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| We investigate the effects of white noise on parametric resonance in X(f> A theory. The 

■ potentiai V(<f>) in this study is ^m 2 (f> 2 + \g4> A + j^(t> 4 ■ An Mathieu-like equation is derived 

and the derived equation is appiied to a partialiy thermaiized system. The magnitudes of 
the amplifications are extracted by solving the equations numerically for various values of 
parameters. It is found that the amplification is suppressed by white noise in almost all the 
cases. However, in some g = cases, the amplification with white noise is slightly stronger 
C^*) ' than that without white noise. In the g = cases, the fields are always amplified. The 

04 ' amplification is maximal at k m 7^ in some g = cases. Contrarily, in the g = 3\/2Am/2 

cases, the fields for some finite modes are suppressed and the amplification is maximal at 
k m ~ when the amplification occurs. It is possible to distinguish by these differences 
^ ' whether the system is on the g = state or not. 

O ' 

cn : 

l . §1. Introduction 

o ■ 

t^J- ' Scalar field theories are widely used to describe many physical systems. One of 

them is A</> 4 theory that is applied to phase transitions. 1 ) These scalar theories are 
also used to describe dynamical phenomena such as nucleation, spinodal decompo- 

"V* ' sition and particle production. It is expected that the results with A</> 4 theory give 

^h. us many indications about the above processes. 

r-j ', Mechanism of particle production is an important topic in various branches of 

physics. The energy transfer from the (false) vacuum to nonzero modes has been 
investigated in terms of particle production. An mechanism of particle production 
is parametric resonance 2 )' 3 ) that the periodic oscillation of an field amplifies other 
fields. Another mechanism is spinodal decomposition that the roll-down of the vac- 
uum from the top of the potential hill with negative mass squared amplifies other 
fields. Particles produced by these mechanisms determine the following evolution 
of the system. An effect by produced particles is the heating of the system as dis- 
cussed in the study of early universe. The similar heating is expected in rebreaking 
processes of the chiral symmetry in high energy heavy ion collisions. 

Parametric resonance in the evolution of the universe 4 ) -24 ) and in chiral phase 
transitions 25 ) -32 ) has been investigated with scalar field theories. An example of an 
oscillator which periodically moves is the condensate of a scalar field. The oscillation 
of the condensate amplifies the other fields, which corresponds to particle production. 
If the expansion of the system is not too rapid, these produced particles will be 
thermaiized by collisions between them. This thermalization corresponds to the 
heating of the universe. In chiral phase transitions, pions (or sigma mesons) will be 
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produced and thermalized. If the system expands rapidly, the produced particles 
may not be thermalized, because the distance between produced particles becomes 
long. Such a case may occur in heavy ion collisions, because the system expands 
one or three dimensionally. 33 ) Therefore, the momentum spectrum may have a peak 
that corresponds to the momentum of the produced particles due to parametric 
resonance. 29 ** 

In a realistic physical phenomenon such as chiral phase transition, friction and 
noise appear because of the existence of an environment in general. For example, the 
linear a model includes only ir and a mesons. However, in a heavy ion collisions, there 
are other particles of different kinds in practice. Then these particles constitute an 
environment. Even when the theory describes the system completely, some fields can 
be regarded as an environment if those fields are projected out. 34 ) -37 ) An candidate 
for an environment in an effective theory at low energies is 'hard modes'; Friction 
and noise appear by integrating the hard modes out if Xcj) 4 theory is regarded as 
an effective theory at low energies and/or if the attention is paid only to the soft 
modes. Therefore, the effects of friction and noise have been investigated in X(p 4 
theory, because the equation of motion has friction and noise terms effectively. 

It has been recognized that noise plays an important role in various fields. An 
attractive example is stochastic resonance 38 - 1 that signal is enhanced by noise. Some 
modes on parametric resonance may be enhanced strongly by noise as found in the 
study of stochastic resonance, though noise disturbs a clear signal in general. The 
effects of noise on parametric resonance have been studied. 21 )- 23 )> 39 ) i n such studies, 
it is found that noise amplifies the fields on parametric resonance. However, it is 
expected that the amplification of a field depends on the form of the coupling term 
to noise. The coupling form between the amplified field <f> and the noise £ is £<^> in 
Ref. 21). However the coupling form in \<j> theory is different. 

In this paper, we study the effects of white noise on parametric resonance in 
A</> 4 theory. We use a classical equation of motion as a practical tool, because 
the amplification of a classical field can be a signal such as particle production 
in quantum mechanical system. Friction and noise terms are attached by hand to 
reveal the role of noise. In sec. |21 the basic equation that includes the effects of 
noise is derived. In subsec. 12.11 the equation of motion for the soft modes is derived 
in the case that the condensate moves periodically near the bottom of the potential. 
In subsec. 12.21 the adopted method in subsec. 12.11 is applied to the A</> 4 theory at 
finite temperature. The derived equation is used to reveal the effects of white noise 
in the subsequent section. In Sec. |HJ the numerical results are presented for various 
values of parameters by solving the equation derived in Sec. 12.21 on the symmetric 
vacuum and the asymmetric vacuum. Sec. |I]is devoted to the conclusions. 
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§2. Formulation 

2.1. Equation of motion with white noise 
The Lagrangian in Xcj) 4 theory is 

C = ^m-^{<p 2 -v 2 ) 2 , (2-1) 

where v 2 can be positive or negative. However, we do not treat the v 2 = case in 
this paper. Another treatment is necessary in such a case, because the mass term 
vanishes. The minimum of the potential is (j> = for v 2 < 0. The potential V(4>) can 
be rewritten as V(<j>) = ^m 2 (j) 2 + \\4> 4 + \Xv 4 , where m 2 = —Xv 2 . On the contrary, 
the minima of the potential are ±v for v 2 > 0. Here, we rewrite the potential 
V{4>) in terms of the new field <fi given by 4> = <j) + v. The potential is rewritten as 
V(4>) = ^m 2 <f> 2 + \g4> i + jA(/> 4 , where m 2 = 2Xv 2 and g = 3Xv. Therefore we use the 
following potential in this paper: 

V{<f>) = \m 2 ^ 2 + \g^> + \x<j>\ (2-2) 

For this system, the Euler-Lagrange equation is 

U(f) + m 2 (j) + g(f) 2 + X(f) 3 = 0, (2-3) 

where □ = d 2 — V 2 . We add friction and noise terms by hand which appear through 
the projection to soft modes, 34 )~ 36 ) the Caldeira-Leggett method, 37 ) the decomposi- 
tion of the effective action 40 ) and so on. Therefore, for slowly varing fields (i.e. soft 
modes) , we use the following equation corresponding to eq. (|2-3p as a practical tool: 

□0 + Ti(t)dt<f> + m 2 + g4> 2 + Xcp 3 = £(t), (2-4) 

where rjdt(f> is the friction term and £ is the noise term. We investigate the systems 
described with eq. (|2-4(l in this paper. 

Because we study the amplification of the finite modes, we divide the field <fi into 
two parts as follows: 

</> = <P + V, (2-5) 
where $ is the spatial average of <j). Substituting eq. ()2-5jl into eq. Q2-4J1 . we obtain 

d 2 <P + r]d t <P + m 2 <P + g<P 2 + A^ 3 

+ Otp + r]d t ip + (m 2 + 2g<P + 3A<2> 2 ) ^ + (g + 3A<2>) ^ 2 + A^ 3 = £. (2-6) 

If ip is sufficiently small at the beginning of the evolution, the lowest order 
equation is 

d 2 t <Z> + r]d t <P + m 2 <P + g<P 2 + A<Z> 3 = f . (2-7) 
In the region that eq. (|2-7j) holds, the equation for ip is 

□V> + v d t ip + (m 2 + 2g<P + 3A<?> 2 ) if) + (g + 3A<2>) ip 2 + A^ 3 = 0. (2-8) 
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We try to solve eq. (|2-8|) with the background field 

The equation of motion for <P near the bottom of the potential with small am- 
plitude is approximately obtained by neglecting <P 2 and ^ 3 terms in eq. (j2-7|) . This 
equation is 

d?& + 7]d t $ + m 2 $ = £. (2-9) 

The solution of eq. (|2-9|) is described as <P = <Pn + ^ih, where is the solution of 
eq. (j2-9j) with £ = and <Pm with £ / 0. $h in the <9 t ?? ~ case is 

^h(*) ~ exp ^— i ^ dsrj(s)\ C + exp (is a + (s)^ + C_ exp cis a_(s)^ 

' ! ' ' ' (2T0) 

where a±(t) = ± (j] 2 /4 — m 2 ) 1//2 and C± are constants. C± are real for r/ 2 /4— m 2 > 
and complex with the relation C+ = C* for ?] 2 /4 — m 2 < 0, because is real. 
^ih is obtained by harmonic analysis. <&m and £ are expanded to Fourier series: 



#ih= ^ ^ne W , (2-lla) 

n=— oo 

oo 

n=— oo 

where w n = 2im/T, n is an integer and T is time interval. Substituting eqs. (|2Tla|) 
and (12-llbj) into eq. (t2~3j) . we obtain 

— u)£ + zryu; n + m z 

if r/ is independent of time. Thus, the approximate solution is expressed with eqs. 
(t2~TT71) and (12- Hal) as follows: 

<f> ~ <2> H + <P IH . (2-13) 

We use this solution to solve eq. (|2-8|) . is evidently a random variable, because 
£ is included in eq. (|2-13|) . In order to obtain the approximate equation for that 
includes the effects of white noise, we replace random variables in eq. Q2-8|) with 
statistical averages. Here, the statistical average is denoted by (•••)• For example, 
in eq. (|2-8|) is replaced with 

Hereafter, we assume that £ is white noise: £ satisfies (£(t)} = and the power 
spectrum 7(u;;£) of £ is a constant I c . The statistical averages (#) and (^ 2 ) are 

(<P) = <2> H , (2T4a) 
<<2> 2 > = (<P h ) 2 + <(<Z>ih) 2 >. (2-14b) 



From eq. (|2-12|) . the relation between the power spectrum / (uj, 'Pm) of <^ih and that 
of £ is 

(uj 2 — m 2 ) 2 + t] 2 uj 2 
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Wiener-Khinchin theorem gives the relation between the power spectrum and the 
correlation function. This theorem holds under the condition that the system is 
stationary, that is, the correlation function depends on only time interval. This 
condition does not always hold in the present case, because (£>h exists. This condition 
holds approximately if the motion of the condensate is sufficiently slow. Moreover, 
we treat the cases that the effects of the initial condition for <^ih are weak. This 
theorem gives the following relation in such a case 



0M*i)«Mt2)> = i c 



duj- 



/-oo (w 2 — m 2 ) 2 + T] 2 U) 2 
We evaluate eq. (|2T6|) for m 2 > -n 2 /4. The integral in eq. 1)2 T6|) is 



(2-16) 



du- 



TT 



oo (a; 2 - m 2 f + r\ 2 u 2 rn 2 r] 
2 v 2 /4) 1/2 \h-t2 



exp 



1 



r?7 |*i - t 2 



x < cos 



[m 



7] sm 



+ 



(m 2 - 77 2 /4) 1/2 |ii-i 2 



2 ( m 2_ ^2/4)1/2 



Therefore, the statistical average of [#ih] is given by 
From eq. (|2T()|) . can be rewritten as 



m 2 r\ 



<1>H = C cos([m 2 — i] 2 exp 



ds rj(s] 



t 



(2-17) 



(2-18) 



(2-19) 



where C and 9 are constants. We obtain 



(2g3> + 3A<2> 2 ) = + -AC 2 exp 

m l r\ 2 



dsr](s] 



to 



+ 2gC exp 



dsrj(s) I cos 



to 



[m 



W4) 1/2 * + 



to 



+ 3AC 2 exp(-/ dsr](s) ) cos 2 (m 2 - ?? 2 /4) 1/2 t + 20 . (2-20) 



Next, we derive the equation for ip when the amplitude of if) is small. Higher 
terms of ip are ignored in such a case. After Fourier transformation, eq. 1)2-8(1 
becomes 

<9 2 V> + r]d t ip + (A? + m 2 + 2 5 # + 3A<? 2 ) ^ = 0, 



(2-21) 

where represents the mode. In order to remove dtip term, we rewrite eq. (|2-21|) in 
terms of the new field Xi which is related to ip as follows: 



ip = exp 



dsr](s) x(f). 



to 



(2-22) 
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We obtain the equation in the dtf]{t) ~ case: 

,2 



dfx + [k +m 



IL + 2g $ + 3\$ 2 ) \ - 0. 



(2-23) 



Replacing 2## + 3A# 2 with (2#0 + 3A$ 2 ) and substituting eq. (|2~2U|) into eq. 
we obtain 

,74) 1/2 



d 2 X + {a(t,k)-2 qi (t) 



cos 



m 



where 



—2q2(t) cos 
a(i, k) = k 2 + m 2 - 



2 (m 2 -77 2 /4) 1/2 t + 2fl]} 



r7 2 3vrA/ c 3 ^ 2 , 

T + + o C Aexp 

4 m^ry 2 



t + e 
x = o, 



dsr](s] 



to 



2q 1 (t) 
252 (*) 



-2pC exp 
-^AC 2 exp 



ds7](s] 



to 



ds7](s] 



to 



(2-24a) 

, (2-24b) 
(2-24c) 
(2-24d) 
r? 2 /4) 1/2 t + 



Rewriting eq. (|2-24a|) in terms of the new variable z given by 2z = (m 2 
9, we obtain 

d 2 z x + [a(t, k) - 2qi(t) cos (2z) - 2g 2 (t) cos (4z) }x = 0, (2-25) 

where a = (4a)/(m 2 -rf/A), q\ = (Aqi) / (m 2 - rf / A) and g 2 = (Aq 2 ) / (m 2 - rj 2 / A) . We 
solve eq. Q2-25|) numerically in the following section. Note that q\ is equal to zero in 
the g = cases. Amplified modes are determined by a for small q± and small q 2 - It 
is found that the amplified modes become soft, because noise increases a. Therefore, 
some modes may vanish because of noise. 

2.2. Application to a Partially Thermalized System 

In this subsection, we treat a partially thermalized system in which the hard 
modes are thermal. This situation may be realized after particle production due to 
spinodal decomposition in early universe and in chiral symmetry rebreaking process 
in high energy nucleus-nucleus collisions. We apply the method described in the 
previous subsection. 

We use A</> 4 theory and divide the field </> into three parts, the condensate <P, soft 
modes (f> s and hard modes 0h. 30 )' 41 ) It is assumed that the hard modes are thermal. 
The thermal average of </> 2 1 n+1 vanishes in the free particle approximation 

(2-26) 



^ +1 ) T = 0, 



where (• • • )t indicates the thermal average to the hard modes and n is a non- negative 
integer. The thermal average of Ac/) 4 potential to the hard modes with eq. ()2-26(l is 



(V(<l>)h 



m 
~~2 



i 2 + |^ + -<* 4 
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= 9(<®t* + \m\T)& + §<2> 3 + \^ + g{^ 8 + \m 2 {T)£ + |^ + ^ 

+ m 2 s $ + g {<p s <P 2 + 4>l$) + ~ ^ 3 s + 6<£ 2 s 2 + 4$tf + 12*&<$) T + <0h>T 

(2-27) 

where m 2 (T) = m 2 + 3A(0 2 1 )t- Kinetic part is expanded in the same way: 

~ (d^<P) T =\d^d^ + \d^ s d^ s + dpfad** + ^(d^ h d^ h ) T . (2-28) 

The (</> 2?1 )t and {d^hd^tphfT terms are independent of <P and <f> s in the massless free 
particle approximation. 42 ) The Euler-Lagrange equation for $ and 4> s after taking 
the thermal average to the hard modes is 

£ + m 2 {T)$ + g$ 2 + A<P 3 + g{<fc) T 

+ + m 2 {T)<j) s + g<j) 2 + A0 3 + 2<?</> s <2> + 3A (<2>0 2 + <2> 2 </> s ) = 0. (2-29) 

Here, the friction and noise terms are added by hand as in the previous subsection. 
The approximate equations for <P and (j) s are, respectively, 

$ + r]$ + m 2 {T)<P + g$ 2 + A<2> 3 + g(<pl)T = £, (2-30a) 

□0 S + r)<j> s + m 2 (T)0 s + (2g<P + 3A<Z> 2 )c/> s = 0. (2-30b) 

In order to remove the g{4>h) r T term in eq. (|2-3Uaj) . we rewrite eq. (|2-3Ua[) in terms 
of the new field E given by = E + <Pq, where is the solution of the equation 
A<£g + g<pQ + m 2 {T)<pQ + ^(0 2 )t = 0. According to the calculation described in 
the previous section, we obtain the equation for <p s by Fourier transformation and 
changing of the variable: 

(2-31a) 
(2-31b) 



(2-31c) 

(2-31d) 
(2-31e) 

(2-31f) 
(2-31g) 



d 2 X + {a T - 2q[ cos(2z) - 2q% cos(4z)} X = 0, 



where 

(j> s (t, k) = exp 



to 



dsrj(s) x(t,k), 



k 2 + M 2 (T)-^ + 



n 2 3vrA/ c 3 ^ 2 . 

+ -C 2 Aexp 



ds r](s 



4 M 2 (T)r? 2 
- -l 

M 2 (T)-rj 2 /A 

qf = -4(g + 3A<2> ) Cexp (-1 ds r,(s^j [M 2 (T) - r, 2 /A\ ^ , 



= —3 AC 2 exp 



2z = M 2 (T) 



f ds r](s)) [M 2 {T)-r, 2 /A\ 1 
to ) 



t + e, 



M 2 {T) = m 2 (T) + 2g<P + 3X^. 
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The remaining task in this subsection is to evaluate m 2 (T) and I c . m 2 (T) includes 
the effects of the hard modes. (0^)t is calculable with the cutoff A between soft and 
hard modes. 

m 2 (T) =m 2 + AT 2 [°° du^—, (2-32) 
2?r Ja/t e u -1 

where the vacuum contribution to (0u)t is dropped. If A is sufficiently small and/or 
T is sufficiently high, eq. 1)2-32(1 is approximately 

m 2 (T) ~ m 2 + ^T 2 . (2-33a) 

On the contrary, if A/T is sufficiently large, 

m 2 (T) ~ m 2 + — IT A exp f-yl/T) . (2-33b) 

I c is related to (£ 2 h )t with eq. (|2-18[1 . We treat the Eb, = case to determine 
the expression of I c . The energy E(S) of 3 for small (j) s is approximately 



£(5) = V 



l -{d t 3f+ l -M 2 (T)3 2 



(2-34) 



where V is the volume of the system. The statistical weight dw is 



M 2 (T)V 

dw = K > exp {-E{3 m )/T) d(d t S m )dS m . (2-35) 

Z7Ti 



Thus, the statistical average ((Sm) 2 ) is 



<( ~ IH ) 2 > = y du; ( ~ IH ) 2 = — (2-36) 
With the help of eq. (|2-18|) . we obtain 

F = ^ (2-37) 

7T|/ 

As shown in Ref. 35), this indicates that the correlation of £ is represented as follows: 

{i{ti)t(t2)) = 2vr/ c <5(t 1 - t 2 ) = ^S(h - t 2 ). (2-38) 
Therefore, the contribution from Sm in a T is r\ independent: 



[M(T)] 2 r/ [M(T)]V 

At high T, eq. (|2-39j) is proportional to (VT)^ 1 and independent of A. On the 
other hand, qj and q 2 go to zero as T goes to infinity. Therefore, the fields are not 
amplified at extremely high T. In a partially thermalized system with white noise, 
eq. (|2-31a|) with eq. (|2-37|) will be solved. 
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§3. Numerical Study 



There are many applications of A</> 4 theory. Inflation and chiral phase transitions 
are interesting examples. In this section, we study the effects of white noise on 
parametric resonance generally in X(j) 4 theory by using scaled parameters. Then, we 
can obtain some indications about above subjects from the results if we want. We 
use the parameters scaled by mass m, for example, k m = k/m. The parameters 
and the variables in the calculations presented below are summarized: Xm = x/ m i 
T m = T/m, r] m = r)/m, g m = g/m, M m (T) = M(T)/m and C m = C/m. Though r] 
or I c cannot be given without clarifying the origin of the dissipation, I c is related to 
r] through the fluctuation-dissipation theorem. Eq. (|2-31aj) is rewritten with these 
parameters for A/T ~ 0: 



9 2 2 Xm + {« T -2gfcos(2z) 



rp 

2q 2 cos 



(4^)}Xr 



0. 



(3-la) 



where 



(M m (T)) 



4 + 4 



km + 



3AT m 3 2 

Ml{T){m3V) + 2 m6XP 



dwr/ m (w) 



WQ 



(M m (T)) 2 -^/4 



T 



g m + 3A 



rn 



C m exp 



dwr] m (w) 



WQ 



(M m (T)) 2 -^/4 



T 
Q2 



3AC^ exp 



WQ 



2z 

2 



dw rj r , 

1/2 : 



(M m (T)) 2 -^/4 



-i 



(Ml{T)-r, 2 jAy IA W + e, 



1 + -\T' m + 2g r , 



m 



+ 3A 



rn 



(3-lb) 



(3-lc) 

(3-ld) 
(3-le) 
(3-lf) 



Here, wq in eq. (|3-lb|) corresponds to zq. It is expected that amplification occurs 



around a T = n 2 if the condition \q[\ 3> \q 2 \ is satisfied, where n is a positive 



, T\ 



integer. Note that a T is equal to or larger than 4 for positive A if (M m (T)Y 



/4 



is positive. In this paper, we treat the cases of (M m (T)) 2 — rf^/A > 0. Eq. (|3-la() is 
solved numerically in the calculations presented below. 

The effective potential is obtained by substituting S = into eq. (|2-27|) . The 
effective potential divided by m 4 with eq. ()2-33a|) is displayed in Fig. ^ for (a) 



g m = 0, A = 10 and (b) g m = 
extremum for g = 3Xv and m 2 



3V2A 



A 



10. 



2 — , 'v — The potential at <P = —v is always 
= 2Xv 2 , which is displayed in Fig. Q(b). Curves 
indicate the effective potential for various values of T m . The critical temperature T c 
that is determined from the effective potential for v 2 > is (2/A) 1 / 2 m. We discuss 
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the amplification in only the cases of T < T c for v 2 > in the calculations presented 
below, because the effective potential is convex anywhere for T >T C . 



In the following subsections, we treat the g m = and g m = 3v ^ 2A cases, because 
g m is exactly zero for v 2 < as illustrated in Fig. ^a) and because g m has a finite 
value for v 2 > as in Fig. ^b). g m is for the Lagrangian defined by eq. (|2T|) , 
In the numerical calculations of time evolution of the soft modes, Xm is set to 1 at 
the initial time and the range of the modes k m is < k rn < 1.5. The condition 

rj <C m is equivalent to < r/ m <C 1. Moreover, we cover two cases: one is that r] m 
is independent of A and the other is that rj m is proportional to A. 34 )' 36 - 1 We set the 
initial amplitude C m to (1 — \ / ^f?>) / \/2\ which corresponds to the inflection point 
of the potential at T m = 0. This initial amplitude is used in subsec. 13.11 and 13.21 
The magnitude of the coupling A is taken to be 0.1, 1 and 10. Then the apparent A 
dependence may be found. The remaining parameter is rrc'V that is related to the 
magnitude of the noise. This parameter is included only in the term coming from 
noise. In this paper, rrfiV is set to 10 if there is no explanation. m 3 V is set to 100 
for comparison. 'Amplification' in the figures displayed below represents the ratio of 
the amplitude for z > 1 to that at the initial time. 

3.1. The g m = cases 

In the case of g m = 0, the q[ term vanishes because of = 0. The equation 
is quite similar to the Mathieu equation, except that a T and are time dependent 
functions. It is determined from the position of (a T ', qj) on the a T — qj plane whether 
the fields are amplified. To compare eq. ()3Taj) with the Mathieu equation, we rewrite 
eq. ()3Ta(l in terms of the new variable z given by 2z = 4z. we obtain 

dl Xm + {a T -2qT C os(2z)} Xm = 0, (3-2) 
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where a T = a T /4 and q\ = q\ '/4. Large amplification is expected if (fi is sufficiently 
large, because a T can be 1. 
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10 20 30 40 50 60 

Time (z) 



Fig. 2. Time evolution of the soft mode \m with k m = 0. The values of the parameters are A = 1, 
r)m = 0.1, T m = 0.1 and C m = (1 - l/^)/^2A. 

Typical time evolution of the soft mode Xm with /c m = is displayed in Fig. [2] 
. The values of the parameters are A = 1, rj m = 0.1 and T m = 0.1. The solution 
for z > 1 is like a sine function, because qj is almost equal to zero when z is 
sufficiently large. Therefore, we can determine the amplitude of Xm for z>l, The 
main purpose of this section is to show the amplification of the soft modes. For 
Z>1, the amplitudes for various values of the parameters can be determined from 
the numerical solutions of time evolution. 

Figure |3] displays the temperature dependence of the amplification with A = 1 
and rj m = 0.1. Each line represents the amplification as a function of k m . It is found 
from this figure that the amplification is suppressed at high T m . This behavior 
is explained by the competition between q\ and a T . On one hand, q\ decreases 
monotonically as T increases; on the other hand, a T does not vary monotonically 
as T increases and goes to 4 as T goes to infinity. Note that the curve of the 
amplification with white noise at T m = is the same curve without white noise at 
T m = 0, because the term coming from white noise vanishes at T m = 0. 

The A dependence is displayed in Fig. 0] with rj m = 0.1 and T m = 0.1. It is 
evident that the fields are amplified for small A and low k m . This dependence is 
most likely explained by the A dependence of a T and qj. a T depends on A explicitly 
and implicitly. Because the A dependence of M m at low T m is weak, the A dependence 
of a T is resemble to the function c\ + C2A, where c\ and C2 are constants. (Because 
C m is (1 — l/v / 3)/v / 2A in this calculation, the term with C m in a T is independent 
of A.) q\ also depends on A weakly through M m . Therefore, the fields are amplified 
for weak A. Moreover, a T increases monotonically as a function of k m . Hence, the 
amplification is strong for small k m . 

Friction dependence is displayed in Fig. |SJ Strong friction make q^ small, 
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because of the exponential factor in . Especially, the peak appears at k m ^ 
for strong friction, while it appears at k m = for weak friction. The suppression 
by friction is strong for the soft modes in this figure. This suppression is probably 
explained by the following fact. The amplification depends not only on the initial 
magnitudes of a T and q%, but also on the time dependences of these parameters, 
while the amplification is determined by the (initial) magnitude of a T and qj in the 
Mathieu equation. 

It seems that white noise always suppresses the amplification. However, it en- 
hances the amplification in some exceptional cases. Figure displays the amplifi- 
cation with A = 0.1, r] m = 0.1 and T m = 0.1. Thin curve shows the amplification 
with white noise and dotted curve without white noise. In the region of k m < 0.1, 
the amplification with white noise is slightly stronger than that without white noise. 
This inversion seems the contradiction to the theory of the Mathieu equation, be- 
cause a T with white noise is always larger than a T without white noise and because 
a T is always larger or equal to 4. In the g = cases, the amplification is determined 
by a T /4 and —qJ/4. Figure [7| displays the time dependence of a T /4 and —q%/4 
with k m = 0, A = 0.1, r\ m = 0.1 and T m = 0.1. It is easily found that — goes 
to zero quickly in a few oscillations of the condensate. This fact implies that the 
magnitude of the amplification is determined not only by the magnitudes of a T and 
qj, but also by the time dependence of them. Then, the inversion does not indicate 
the contradiction and happens in some cases. 

Figure |H] displays the amplification for m 3 V = 10 and 100 with A = 1 and 
T m = 0.1. Figure |H1 (a) displays the amplifications with r] m = 0.1. Figure |H1 (b) 
displays with r\ m = 0.2. As shown in Fig. El (a), the amplification for rrc'V = 100 
is slightly stronger in wide range of k m than that for m 3 V = 10, because large 
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Fig. 5. Friction dependence of the amplifica- 
tion with A = 1 and T m =0.1. The square 
designates data for <r) m = 0.05, the circle 
for r)m =0.1, the triangle for rj m = 0.2 and 
the diamond for r\ m = 0.3. 



< 




with white hoise 
without white noise 



km 



Fig. 6. Comparison between the results with 
white noise and those without white noise. 
The values of the parameters are A = 0.1 
and T m = 0.1. The square designates 
data for 'with white noise' , the circle for 
'without white noise'. In the region of low 
km, the amplification with white noise is 
stronger than that without white noise. 
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Fig. 7. Time evolution of a T /4 and -gJ/4 with k m = 0, A = 0.1, Tj m = 0.1 and T m = 0.1. a T /4 
and — q[ /4 decrease quickly. — <jg /4 approaches to zero asymptotically. 



volume reduces the magnitude of the white noise. However, the magnitude of the 
amplification for m 3 V = 100 is smaller than that for m?V = 10 in the region of small 
k m as shown in Fig. El(b). As discussed in the previous paragraph, this inversion is 
probably caused by the quick decrease of ■ 

In any case, the fields are amplified in the soft modes. The amplification de- 
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creases monotonically as a function of k m in the region of large k m and goes to 1 
asymptotically as k m goes to infinity 



1 .4 I 1 1 1 1 1 1 1 r 




j I I I I I L 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
X 

Fig. 9. Amplification in the cases that r\ m is proportional to A. The values of the parameters are 
km = and T m = 0.1. The friction is parametrized as r\ m — T] m oX. The square designates data 
for r)mO = 0.1 and the circle for 7/ m o = 0.2. 



rj m in the above calculations is independent of A. In the last part of this subsec- 
tion, we investigate the amplification in the A dependent cases that r\ m is proportional 
to A as r\ m = r/ m oA, where rj m o is independent of A. The A dependence of the ampli- 
fication for rj m Q = 0.1 and 0.2 are displayed in Fig. |UJ In Fig. |5J the amplification is 
strong in the weak coupling region and weakens quickly with the increase of A. It is 
clear that strong friction causes this suppression, because the friction decreases q£ . 
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3.2. The g m = cases 

In this subsection, we investigate the amplification in the g m = cases. q\ 
is not zero in the present cases. [See eq. (j3-lcj) .] The ratio (qiV^JI is larger than 
1 in almost the time, especially, in the late time. Therefore, q[ plays an important 
role on the amplifications of the fields . 




Figure ^3 displays the temperature dependence of the amplification with rj m = 
0.2 for (a) A = 1 and (b) A = 10. The amplification in the region of low k m is 
suppressed and this suppression is strong for large A. Therefore, the amplification 
may occur only for weak A at finite temperature. Dips are found around k m = 0.6 
at T m = in the present cases, while there is no dip in the g m = cases. Note that 
there is no effect of white noise at T m = 0. The curves accord well in the region of 
k m > 1 for A = 1 and tend to 1 as k m increases. This tendency in the region of large 
k m is also displayed for A = 10. This can be explained by the k m dependence of a T 
[eq. (|3-lb|) ]. We remember that the Mathieu equation has the following form: 

d 2 z y{z) + {a-2q cos(2z)) y(z) = 0. (3-3) 

The stable region of the Mathieu equation spreads according as a increases. There- 
fore, the solution of the Mathieu equation is similar to a sine function for large a, 
because the cosine term in the Mathieu equation is negligible. The amplitude of the 
field is equal to the initial value 1, because no amplification occurs for large a in the 
Mathieu equation. The time dependence of the solution of eq. (|3-la|) for large k m is 
explained on similar lines, because eq. Q3TaJI is similar to eq. ()3-3|) . In Fig. El the 
amplification in the region of low k m for T m ^ is weaker than that for T m = 0. This 
behavior is a nontrivial result from the equation similar to the Mathieu equation. 
The amplification for A = 10 does not decrease monotonically as the temperature 
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Fig. 11. Coupling dependence of the amplification with r/m = 0.2. Figure till (a) is depicted for 
T m = 0, Fig. EI](b) for T m = 0.1 and Fig. HTICcl for T m = 0.2. The square designates data for 
A = 0.1, the circle for A = 1 and the triangle for A = 10. 



increases. This reduction is most likely explained by the oscillation of the fields and 
the time dependence of the parameters, a T , q[ and qj. 



3.5 I 1 1 1 1 1 1 1 1 3.5 




Fig. 12. Comparison between (a) the coupling dependence of the amplification without white noise 
and (b) that with white noise. The values of the parameters are T m — 0.2 and r} m =0.1. The 
square designates data for A = 0.1, the circle for A = 1 and the triangle for A = 10. 



The amplification for various values of A is displayed in Fig. ^Jwith rj m = 0.2 
and various values of T m . The amplification at T m = is displayed in Fig. ITTT a). 
The curves in this figure are the same, because eq. (|3Tbj) is A independent at T m = 0. 
Indeed, AC^ is A independent, because the initial amplitude of the condensate C m 
is proportional to A -1 / 2 . Moreover, the coupling g and the mass squared m 2 are 3Xv 
and 2Xv 2 , respectively. Thus, is A independent, because <Po satisfies the equation 
$o + 3t>$o + 2v 2 <P = 0. Therefore, a T , qj and q\ are A independent at T m = 0. 
We also remember that rj m is A independent. (This value is given by hand in the 
present numerical calculations.) From this consideration, we find that eq. (|3-la|) at 
T m = is independent of A. (This coincidence is artificial in the sense that the result 
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depends on the initial amplitude of the condensate and the A dependence of friction.) 
The A dependences at T m = 0.1 and T m = 0.2 are displayed in Fig. fTIT b) and Fig. 
Illf c) respectively. The amplification weakens in the region of low k m according as 
A increases. Especially, the amplification is strongly suppressed for A = 10. 

Figure PHZl displays the amplification for various values of A with r/ m = 0.1 and 
T m = 0.2. Figure lTW a) displays the amplification without white noise and Fig. I12f b) 
with white noise. The strength of the amplification increases with the coupling A in 
Fig. 112( a). while that decreases with the coupling A in Fig. I12f b). The A dependence 
of the amplification is opposite. This result implies that white noise can change the 
growth of the fields. (Therefore the particle distribution is probably different.) 




Figure IT31 is depicted to present the time dependence of a T without white noise 
for various values of A with 7] m = 0.1 and T m = 0.2. a T for A = 10 is largest in 
the early time, but it is smallest in the late time. (Note that the curves in Fig. 
intersect at different points.) The z dependence of the parameters implies that the 
field can grow even for large A as shown in Fig. I12f a), 

Figure ITU displays the amplification for some r] m with A = 0.1 and T m = 0.1. 
The strength of the amplification decreases with the friction rj m . Apparently, the 
friction makes the amplification weak. The reason is the same as discussed in the 
previous subsection. 

Figure ED displays the amplification for different volumes with (a) rj m = 0.1 
and (b) r\ m = 0.2. The parameters are A = 1 and T m = 0.1. In both cases, the 
amplification for m 3 V = 10 is weaker than that for m 3 V = 100. (The amplification 
for m 3 V = 10 is slightly stronger than that for m 3 V = 100 around k m = 0.7.) a T 
for rr?V = 100 is smaller than that for m 3 V = 10, because the magnitude of white 
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noise is proportional to V as in eq. (|2-38|) . Therefore, the fields are amplified 
strongly for large V . 
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Fig. 16. Amplification in the cases that r) m is proportional to A. The values of the parameters are 
km = and T m — 0.1. The friction is parametrized as rjm = ??moA. The square designates data 
for rjmo = 0.1 and the circle for rj m o = 0.2. 

Finally, as shown in Fig. 1141 the amplification for various values of A is displayed 
in Fig. ED with k m = and T m = 0.1 in the case that the magnitude of the friction 
is proportional to A. The friction is parametrized as r\ m = r] m Q\ as in Sec. 13.11 while 
the friction is A independent in the previous calculations in this subsection. The A 
dependence of the amplification in the g m = 3\/2A/2 cases is similar to that in the 
g m = cases. The amplification depends on the magnitude of the friction strongly, 
because the friction affects the coefficients, a T , q[ and through the exponential 
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factors. Therefore, the strength of the amplification decreases with A as shown in 
Fig. d 

§4. Conclusions 

We investigated the effects of white noise on parametric resonance in A</> 4 theory. 
The potential V{4>) in this study is ^m 2 2 + ^g<fi 3 + |A^ 4 . The effects of white noise 
for the fluctuation field ip defined in eq. (|2-5f) are included through the motion of 
the condensate. White noise modifies the mass term in the equation of motion for 
the soft modes. 

It was found from the present numerical study that the amplitudes of the fields 
for the soft modes are amplified due to parametric resonance in both the g = and 
g = 3 X^ m cases. The white noise suppresses the amplification in almost all the 
cases. However, in some cases, the amplification with white noise is slightly stronger 
than that without white noise. The strength of the amplification decreases with the 
temperature, the friction and the coupling A. This result lead to an opposite conclu- 
sion: white noise suppresses the amplification in almost all the cases in A</> 4 theory, 
while white noise enhances in some other theories. The suppression is caused by the 
<P 2 (j) s term in the equation of motion for the soft modes in A0 4 theory. Therefore, 
the form of the coupling between the condensate and finite modes is essential to the 
amplification due to parametric resonance with noise. 

In the present numerical study, we found two different points between the results 
in the g = cases and those in the g = 3 X^ m cases: 1) The fields are always 
amplified in the g = cases. Contrarily, the fields are suppressed around k m = 0.6 
in the g = 3v ^ m cases. Thus, the relative ratio of the amplitudes for different 
modes is emphasized in the latter cases. In other words, the momentum dependence 
of the field amplitudes is emphasized. 2) The amplification of the field is maximal 
at k m ^ in some g = cases, while that is maximal at k m ~ in g = 3v ^ m 
cases when the amplification occurs. It is possible to distinguish by these differences 
whether the system is on the g = state or not. 

In this paper, the effects of the linear term of £ are not taken into account, 
because the noise is averaged. Back reaction and effects of the nonlinear terms are 
also ignored. We would like to include these effects in the future studies. 
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